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Abstract 

We apply the causal interpretation of quantum mechanics to homogeneous and isotropic quantum cos- 
mology where the sources of the gravitational field are either dust or radiation perfect fluids. We find 
non-singular quantum trajectories which tends to the classical one when the scale factor becomes much 
larger then the Planck length. In this situation, the quantum potential becomes negligible. There are no 
particle horizons. As radiation is a good approximation for the matter content of the early universe, this 
result suggests that the universe can be eternal due to quantum effects. 
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1 Introduction 



The appearance of initial singularities in the classical cosmological models which better describe the universe 
we live in constitutes a big puzzle to all cosmologists. Until now, singularities are out of the scope of any 
physical theory. If we assume that a physical theory can describe the whole Universe at every instant, even 
at its possible moment of creation (which is the best attitude because it is the only way to seek the limits of 
physical science), then these classical singular points must be avoided. Indeed, no one expects that classical 
general relativity continues to be valid under extreme situations of very high energy density and curvature. 
In particular, it is very plausible that quantum gravitational effects become important under these conditions, 
eliminating the singularities that appear classically. To see if this is indeed the case, we should construct 
a theory of quantum cosmology. However, any quantum theory when applied to cosmology presents new 
profound conceptual problems. How can we apply the standard probabilistic Copenhaguen interpretation to 
a single system as the Universe? Where in a quantum Universe can we find a classical domain where we 
could construct our classical measuring apparatus to test and give sense to the quantum theory? Who are 
the observers of the whole Universe? This is not a problem of quantum gravity alone, because there is no 
problem with the concept of an ensemble of black holes and a classical domain outside it, but it is specific 
of quantum cosmology. As we cannot apply the Copenhaguen interpretation to quantum cosmology, we will 
adopt an alternative non-probabilistic interpretation, which circumvents the measurement problem because it 
is an ontological interpretation of quantum mechanics: it is not necessary to have a measuring apparatus or a 
classical domain in order to recover physical reality. It is the causal interpretation of quantum mechanics |) . 

In this letter, we apply the causal interpretation of quantum mechanics to some specific quantum states 
already found in the literature |3|, and show that the Bohmian (quantum) trajectories are non-singular, 
tending to the classical trajectories when the scale factor is large. The quantum potential is the responsible 
for this behavior. When the scale factor is small, the quantum potential becomes large creating an effective 
repulsive force avoiding the singularity. When the scale factor is large, the quantum potential becomes negligible 
and the classical potential dominates. The minisuperspace models we study are constituted of a Friedman- 
Robertson- Walker (FRW) metric with either dust or radiation perfect fluids ||, Q]. The ADM quantization 
procedure is performed because in these cases there is a preferable time variable which rends the quantum 
equations into a very simple Shroedinger form. The quantum solutions obtained are gaussians §§. The dust 
case is rather academic but the radiation case is important because it can mimic quite well the matter content 
at the very early universe. In all cases we obtain non-singular eternal models without any particle horizon. 

This letter is organized as follows. In the next section we make a summary of the causal interpretation 
and its application to quantum cosmology. In section 3 we abridge the results of Refs. which concern 

this letter. In section 4 we apply the causal interpretation to the solutions presented in section 3, obtaining 
the Bohmian trajectories. We compare our results with those of Refs. || ij. We end with comments and 
conclusions. 

2 The causal interpretation of quantum mechanics 

In this section, we will review the causal interpretation of quantum mechanics. Let us begin with the 
Schrodinger equation, in the coordinate representation, for a non-relativistic particle with the hamiltonian 
H =p 2 /2m + V{x): 



ift— ^ = h^-V 2 + V(x)Mx, t). (1) 



-2 

dt 1 2m 

Writing \& = Rexp(iS/h), and substituting it into ([!]), we obtain the following equations: 

dS (V5) 2 H 2 V 2 R 



dt 2m 2m R 

dR 2 _ , „, VS 



V = 0, (2) 



+ V.(i? 2 — ) = 0. (3) 
at m 

The usual probabilistic interpretation takes equation (||) and understands it as a continuity equation for 
the probability density R 2 for finding the particle at position x and time t if a measurement is performed. 
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All physical information about the system is contained in R 2 , and the total phase S of the wave function is 
completely irrelevant. In this interpretation, nothing is said about S and its evolution equation (||). However, 
examining equation (^|) , we can see that VS 1 / m may be interpreted as a velocity field, suggesting the identifi- 
cation p = VS 1 . Hence, we can look to equation (|^) as a Hamilton- Jacobi equation for the particle with the 
extra potential term — h 2 V 2 R/2mR. 

After this preliminary, let us introduce the causal interpretation of quantum mechanics, which is based on 
the two equations (Q) and (^), and not only in the last one as is the Copenhaguen interpretation: 

i) A quantum system is composed of a particle and a field 'J (obeying the Schrodinger equation (jl])), each 
one having its own physical reality. 

ii) The quantum particles follow trajectories x(t), independent on observations. Hence, in this interpretation, 
we can talk about trajectories of quantum particles, contrary to the Copenhaguen interpretation, where only 
positions at one instant of time have a physical meaning. 

iii) The momentum of the particle is p — VS. 

iv) For a statistical ensemble of particles in the same quantum field W the probability density is P = R 2 . 
Equation (|J) guarantees the conservation of P. 

Let us make some comments: 

a) Equation (||) can now be interpreted as a Hamilton- Jacobi type equation for a particle submited to an 
external potential which is the classical potential plus a new quantum potential 

Q = (4) 
^ 2m R K 1 

Hence, the particle trajectory x(t) satisfies the equation of motion 

m^f=-W-VQ. (5) 

b) Even in the regions where \& is very small, the quantum potential can be very high, as we can see from 
equation (|4|). It depends only on the form of not on its absolute value. This fact brings home the non-local 
and contextual character of the quantum potential^ This is very important because Bell's inequalities together 
with Aspect's experiments show that, in general, a quantum theory must be either non-local or non-ontological. 
As Bohm's interpretation is ontological, it must be non-local, as it is. The quantum potential is responsible 
for the quantum effects. 

c) This interpretation can be applied to a single particle. In this case, equation (||) is just an equation to 
determine the function R, which forms the quantum potential acting on the particle via equation (^J) . It is not 
necessary to interpret R 2 as a probability density and hence it may not be normalizable. The interpretation 
of R 2 as a probability density is appropriate only in the case mentioned in item (iv) above. The causal 
interpretation is not, in essence, a probabilistic interpretation. 

d) The classical limit is very simple: we only have to find the conditions for having Q — 0. 

e) There is no need to have a classical domain because this interpretation is ontological. The question on 
why in a real measurement we do not see superpositions of the pointer apparatus is answered by noting that, in 
a measurement, the wave function is a superposition of non-overlaping wave functions. The particle will enter 
in one region, and it will be influenced by the unique quantum potential obtained from the sole non-zero wave 
function defined on this region. The particle cannot jump to other branchs because it cannot pass through 
nodal points of the wave function. 

In the next section we will perform an ADM quantization of Friedman-Robertson- Walker minisuperspace 
models containing either dust or radiation perfect fluids. In these cases, a prefered time variable can be selected 
in terms of potentials of the velocity field of the fluids. As a time variable will be chosen before quantization, 
the quantum equations will be exactly like the Schrodinger equation. Instead of the particle position, as in the 
above example, the single degree of freedom of these quantum models will be the scale factor of the universe. 
The interpretation of the quantum solutions will run analogously to what was described above. The time 
evolution of the scale factor will be different from the classical one due to the presence of an extra quantum 
potential term in the modified Hamilton- Jacobi equation it satisfies. 

1 This fact becomes evident when we generalize the causal interpretation to a many particles system. 
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3 The ADM quantization of dust and radiation minisuperspace 
models 

In this section we present the minisuperspace models which we will analyze using the causal interpretation. 
These models are obtained from the quantization of the dust and radiation fluids through the ADM prescription. 
For the quantization of the radiation filled FRW model we follow Ref . Q . We start with the line element 

ds 2 = -N(t) 2 dt 2 + a 2 (t)a lJ dx l dx 1 . 

where is the metric of constant curvature three-surfaces. The full action for a perfect fluid is given by 

S= I d A xs/^g(R w + p) + 2 / d 3 xVhh i:i K ij , (6) 

JM JdM 



where p is the pressure, R^ is the scalar curvature, h%j is the 3-metric on dM, K % i is the second fundamental 
form of the boundary, and we choose c = I&ttG = 1. 
The action (ffl) can be reduced to 



S r = J dt ap a - tpp v + Sps - NH 

where we are using Schutz's fluid variables ip, A, 7, fi, S |§, in terms of which the four- velocity of the fluid is 
written as 

U v = -{d l/ tp + Xd vl + Od u S). (7) 
A* 

The quantity fi is the specific enthalpy. The super-Hamiltonian TL takes the form 



H 



= -(^ + 6ka] +p i / 3 a~ 3 e s . 



\24a 



For details, see Refs. §, f|. 

Performing an ADM reduction using the conformal-time gauge N = R in the radiation case, the reduced 
action becomes 



dt 



a Pa 



24 



-I- 6ka' 



where k = +1, or — 1, for spherical, flat, or hyperbolic spacelike sections of the three space with metric er^ , 
respectively. The Hamiltonian in the reduced phase space, which has only a as a degree of freedom, takes the 
very simple form 

n 2 

H = t^ + Qka 2 . 
24 

The classical solutions are: 

!sin£, for k = 1 

t, for k = 

sinhi, for k = — 1 

which are the well known solutions for radiation in conformal time. 
The quantized Hamiltonian, in units where K = 1, is 



(8) 



H 



1 d 2 



6ka 2 



24 da 2 

As a > 0, the requirement that H must be self-adjoint leads to the restriction 

V>'(0) = oeip(0) 



(9) 



(10) 



on the wave function ip, where a is a parameter in the interval (—00,00]. Two solutions will be obtained, one 
for a = and the other for a = 00. 
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The propagators for the Hamiltonian (^) are well known. However, we have the extra constraint that 
a > 0. For this reason, the Hilbert space is restricted to functions in L 2 (0,oo). For the case where a = the 
propagator is 

( a , a', t) = G(a, a',t) + G(a, -a', t), (11) 

where G(a,a',t) is the usual harmonic oscillator propagator for a system with mass m = 12 and angular 
frequency u> = \fk. Seting a = in Eq. (O), we take a gaussian wavepacket as the initial normalized state, 



1/4 

— I exp(-/3o 2 ), 



(12) 



where (3 = b + iB, B and b are real, and b > 0. The wave function at time t given by Eqs. (O) and (03) is 



1/4 



6v^ 



n 1/2 



i(y/kt)[0tim(Vkt) - 6iVk] 



(13) 



x exp ■ 



tan(VA:t) 



(Vfct)[/3tan(-\/fct) - 6«VI] 



The expectation value for the scale factor a can be computed from dT3|), and is 



1 / 2 

12 V TTfe 



6 2 sin 2 t + (6 - B tan t) 2 cos 2 1 



^6 2 ^ 2 + (6 -~Bffi 



6 2 sinh 2 t + (6 - B tanh i) 2 cosh 2 i 



for k = 1 
for fc = 
for fc = -1 



For a = oo the propagator is 

G (//) (a, a', t) = G(a, a', t) - G{a, -a', t). 
We take as the initial state the wave packet 



4 n) (R)=[-\ Rexp(-/3R'). 



8b 



1/4 



The evolution of ( |16|) governed by the propagator (|15|) is 

xi?exp 



-i -3/2 



6i-\/fc 



tan(Vfci) 

6i\/fc 



(14) 



(15) 



(16) 



(17) 



ta,n(Vkt) \ cos 2 {^f~kt)[f3 ta,n{\fkt) — 6iVk] 
The expectation value for the scale factor a with the wavefunction ( |l7j ) is 

(a)^ = 2<a)W. 



i? 2 



(18) 



We now turn our attention to the dust filled minisuperspace model presented by Gotay and Demaret 
Once again Schutz's variables and a FRW metric are used. In this case the field <p in Eq. (Q) is the only velocity 
potential which is non null. The time to be chosen will be the proper time of the dust particles defined by 
ip = —t, which is equivalent to choose N = 1. Seting p — and fi = 1, the ADM reduced super-Hamiltonian 
becomes (see Ref. ||) 



n 



A. 

24a 



+ 6ka - p^ 



The reduced Hamiltonian, with the above choice of time, is 

P 2 

H(a,p a ) = ^ + 6ka. 
24a 
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If we perform the canonical transformation 



x = ^V6a 3 / 2 , 



\/6 _! 



the Hamiltonian takes the simple form 



H(x,p x ) = V \ + Kx 2 /\ 



where K = |6 3 / 2 fc. Taking k = 0, the classical solution is simple to obtain. It is x{t) — t or equivalently 
a(t) oc i 2 / 3 . The quantized Hamiltonian for k — is, in units where h = 1, 

d 2 



H 



dx 2 



Once again, the requirement of sclf-adjointness of H yields the boundary conditions ip'(0) — aip(0). Choosing 
a = 0, we can evolve the initial gaussian wave function 



■ 



exp(— (3x 2 ), 



(19) 



by applying the corresponding propagator, analogously to what was done in the radiation case. As before, 
ft = b + iB with b > 0. The wave function at time t is given by: 



il>{x,t) 



1/4 



(1 + 4if3t)~ 1/2 exp 



-px 2 



1 + Ai(3t 



The expectation value for x can easily be computed, and is 



We will now turn our attention to the causal interpretation of the wavefunctions V^^i V* i an< ^ 



(20) 



(21) 



4 The causal interpretation 

The causal interpretation of the above minisuperspace models is straightforward. Let us begin by the dust 
field. In this case, the quantum trajectories are solutions of the following differential equation: 



P., 



1 

r 



dSd 
dx ' 



where Sd is the phase of the wave function (|20|). The general solution is: 

x(t) = x [I6b 2 t 2 + (1 - ABtf] 1/2 , 



(22) 



(23) 



where x$ is an arbitrary positive integration constant. That the mean value of x founded in the previous 
section is, apart from the integration constant, the same function of time as the solution ( pi| ) is not surprising. 
Mean values in the causal interpretation are the same as in the usual interpretation if we also assume that 
the amplitude squared of the wave function, R 2 {x,Q), represents a probability density distribution of initial 
conditions of the quantum trajectories. Then, we obtain from Eq. (fL9h 



(x(t)) = 



1/2 



-26a; 



2 °x [l6b 2 t 2 + (1 - ABt) 2 ]^ 2 dx 



= (27rb)- 1/2 [16b 2 t 2 + (l-4Bi) 2 ] 1 / 2 
which is the result (p|) of the previous section. 



(24) 
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It can be seen from Eq. (|23|) that no quantum trajectory is singular. The scale factor is never zero for any Xq 
greater than zero. Also, the trajectories approach the classical one for large \t\ (remember that x(t) oc a 3 ^ 2 (t)). 
This behavior can be seen by examining the quantum potential, 

Qd ~ 2b [{l-ABt) 2 + m 2 t 2 ] 2 ■ ^ 
Along the quantum trajectory (^3|), the quantum potential turns out to be: 

_ 2bx 2 (l - 2bx 2 ) 

Qd{t) xHt) ' (26) 

where x(t) is given by (p3|). 

From the above equation we can see that the quantum potential goes to zero when \t\ becomes large but it 
is positive of order b when \t\ is small^. In this situation, it works like a repulsive force around the region x = 0. 
Figure 1 shows the curves x(t) and Q(t) (along the trajectory) versus t. They represent classical universes 
contracting from infinity to a minimum size, where its behavior is not classical, and then expanding to infinity, 
getting classical again as the scale factor becomes large. These models have no particle horizon. The integral 
/!«, a _1 (i')di' oc £ X - 2 / 3 {t')dt' diverges. 




Fig. 1: A typical trajectory for x(t) and its corresponding quantum potential (dashed line). We see that when 
x approaches the classical singularity the quantum potential becomes large. 

Let us now turn our attention to the more realistic case of radiation. Again, the application of the causal 
interpretation is straightforward. We have to take the phases of the wave functions ( |l^ , |l7| ) for the cases 
k = 0, k ± 1, calculate their derivatives with respect to a, and equate the result with p a = 12a. The solutions 
of these first order differential equations are: 

{yV sin 2 1 + (6 - B tant) 2 cos 2 t for fc = 1 

^bH 2 + (6 - Btf for fc = . (27) 

b 2 sinh 2 1 + (6 - B tanh i) 2 cosh 2 t for fc = -1 

2 There is nothing special with the choice xq = \j\p2h. The quantum potential is zero along the trajectory but the quantum 
force F d {t) = —d x Qd is not. 



G 



where ao is a positive integration constant. Like before, they have the same functional behaviour as the mean 
values encountered in the previous section for the reasons already explained in the dust case. None of these 
solutions reach the singular point a = .The quantum potentials are given by 



Q = 'ib 



[fc(bsin 2 (t)-72a 2 ) + (6~B t.an(t)) 2 cos 2 (t)] 
[(6-Btan(i)) 2 cos 2 (i)+6 2 sin 2 (t)] 2 

[6(frt 2 -72q 2 ) + (6-Bt) 2 ] 
[(6-Bt) 2 +6 2 t 2 ] 2 

[6(6 sinh 2 (t) -72a 2 ) + (6-B tanh(t)) 2 cosh 2 (t)\ 
[(6-Btanh(t)) 2 cosh 2 (t)+b 2 sinh 2 (t)] 2 



for k = 1 
for k = 
for k = — 1 



At the trajectories, the quantum potentials for k = 0, 1, —1 are the same, and equal to 

36ag(l - 12bal) 



Q 



(28) 



(29) 



where a(t) is given by (27). 

For k = and fe = — 1 the models are qualitatively similar to the dust case. It can be seen from Eq. ( p?| ) 
that these universes contract classically from infinity to a minimum size, where their behaviors are not classical, 
and then expand to infinity, getting classical again as the scale factor becomes large. The case k — 1 deserves 
special attention. Examining the case with B = and b > 6 (the general case is qualitatively similar), we can 
see that the quantum trajectory oscillates between the minimum value 6ao and the maximum value a$b. Their 
ratio is b/6. The ratio between the quantum and classical forces at these points are — & 2 /36 for the minimum 
and — 36/fo 2 for the maximum. The universe we live in is large and classical, and it must have gone through 
a contracted phase where nucleosynthesis took place. This means that the ratio between the maximum and 
minimum values of the radius of the universe (which is the inverse ratio of the respective temperatures at these 
epochs since the fluid is radiation) must be at least of the order of 10 10 . Hence b must be greater than 10 10 , 
yielding a very flat initial gaussian wave function. This ensures that the universe had a nucleosynthesis era 
with a classical behaviour since then. 

These models do not have particle horizons. In this case the particle horizon is proportional to the integral 
J_ a~ 1 (t')N(t')dt' which in the gauge N = a we are using evidently diverges. 

Figures 2 and 3 show the scale factor and the quantum potential for k = and k = — 1. When a is large, the 
quantum potentials go to zero, while when a approaches its minimum size, they become important, creating an 
effective repulsive quantum force around this region. Figures 4 and 5 show the scale factor and the quantum 
potential for k = 1 and for different values of b. Note that for larger 6's the amplitude of oscillation also 
becomes larger. 




40 60 



Fig. 2: Typical quantum trajectory with its corresponding quantum potential (in dashed lines) for the 

radiation fluid model for the case where k = 0. 
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Fig. 3: Typical quantum trajectory with its corresponding quantum potential (in dashed lines) for the 

radiation fluid model for the case where k = — 1 . 




Fig. 4: Quantum trajectories for different values ofb and the same initial condition a(0) for the radiation 
fluid model and k = 1 . The larger the value of b the larger is the amplitude of oscillation. 
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Fig. 5: Quantum potentials computed along the quantum trajectories shown in Fig. 4. The dashed quantum 
potential corresponds to the dashed trajectory in Fig. 4. We can see that the value of the quantum potential 
increases as we approach a classical singularity, being higher when b is higher. 

5 Conclusion 

In this letter we have obtained a class of non-singular cosmological models without particle horizons by applying 
the causal interpretation to some quantum states of the universe already obtained in the literature ^ Q and 
presented in section 3. We have shown that the Bohmian trajectories are the same functions of time as 
their corresponding mean- values obtained in Refs. H which makes use of the conventional Copenhaguen 
interpretation. However, it seems to us that only in the causal approach we can arrive at definite conclusions 
about the existence of singularities. It is not because the mean value of a variable is different from zero at all 
times that this variable cannot be zero sometime. Nothing forbids that this be the case for the single universe 
we live in. The causal interpretation applied to this problem states that each individual trajectory of these 
quantum states is not singular, a much more stronger and valuable result. This means that if the universe is 
in one of these quantum states, then it is indeed non-singular because there is no quantum trajectory with 
singularities. 

For the more realistic radiation fluid, we have obtained singularity- free models, without particle horizons, 
which approach the classical behavior as the universe expands. For flat and negative curvature three-surfaces 
this is the situation in all cases. The universe contracts classically from infinity to a minimum size, where its 
behavior is not classical, and then expands to infinity, getting classical again as the scale factor becomes large. 
For the positive curvature case, the classical behavior is achieved only for some values of the parameters. In 
this case, we would have an eternal periodic universe which is classical when it is large. 

These examples show that quantum gravitational effects can indeed prevent the formation of cosmological 
singularities but this result can only be stated with strength along the lines of the causal interpretation. It 
should be interesting to investigate if the results presented here are stable under small perturbations. Then we 
will have to face new technical and interpretational problems. But this is another story. 
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